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In this article, we perforin a systematic study of the mass spectrum of the axial- 
vector hidden charmed and hidden bottom tetraquark states using the QCD sum rules, 
and identify the Z + (4430) as an axial- vector tetraquark state tentatively. 
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' p ,■ 1 Introduction 

The Babar, Belle, CLEO, DO, CDF and FOCUS collaborations have discovered (or con- 
■ firmed) a large number of charmonium-like states, such as X(3940), X(3872), Y(4260), 

y(4008), ^(3940), V(4325), V(4360), y(4660), etc, and revitalized the interest in the 
spectroscopy of the charmonium states pjj [2j [3j HI O Many possible assignments for 
those states have been suggested, such as multiquark states (irrespective of the molecule 
If} [ type and the diquark-antidiquark type), hybrid states, etc [U 121 El [U [5] • 

The Z + (4430) observed in the decay mode ip'ir^ (B — > iI)'-k + K) by the Belle collab- 
oration is the most interesting subject 0[8]. We can distinguish the multiquark states 
from the hybrids or charmonia with the criterion of non-zero charge. The Z + (4430) can't 
be a pure cc state due to the positive charge, and may be a ccud tetraquark state. How- 
ever, the Babar collaboration did not confirm this resonance [9j. Furthermore, the two 
resonance-like structures Z(4050) and Z(4250) in the tt + Xc1 invariant mass distribution 
near 4.1 GeV are also particularly interesting [10] . Their quark contents must be some 
special combinations of the ccud, just like the Z + (4430), they can't be the conventional 
mesons. There have been several theoretical interpretations for the Z + (4430), such as 
the hadro-charmonium resonance [IT] , the S-wave threshold effect [12] , the molecular 
D*D 1 (D' 1 ) state [I3[ia[ia[IS[I71[l8l[Il[20],the tetraquark state [2JJ [22l [Ml IM 123 [26] , 
the cusp in the D*D\ channel [27], the radially excited state of the D s [28J, the pseudo- 
resonance structure [29], etc. 

In Refs.|30[ 131]. we assume that the hidden charmed mesons Z(4050) and Z(4250) are 
vector (and scalar) tetraquark states, and study their masses using the QCD sum rules. 
The numerical results indicate that the mass of the vector hidden charmed tetraquark state 
is about Mz = (5.12 ±0.15) GeV or (5.16 ±0.16) GeV, while the mass of the scalar hidden 
charmed tetraquark state is about Mz = (4.36 ± 0.18) GeV. The resonance-like structure 
Z(4250) observed by the Belle collaboration in the exclusive decays B° — > K~ir + Xci can 



1 E-mail, wangzgyiti@yahoo.com.cn. 

2 There have been many theoretical works on the X, Y, Z hadrons, it is difficult to cite all of them, 
we prefer the comprehensive review articles [T] [21 [3l [4] [5], where one can find the original literatures. 
On the other hand, one can consult Ref. 6 for a concise review of the experimental situation of the new 
charmonium-like states. 
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be tentatively identified as the scalar tetraquark state [31]. In Refs.[32| 133]. we study the 
mass spectrum of the scalar hidden charmed and hidden bottom tetraquark states in a 
systematic way using the QCD sum rules. In Ref.[34J, we study the mass spectrum of the 
vector hidden charmed and hidden bottom tetraquark states systematically. Recently, the 
hidden charmed and hidden bottom tetraquark states are studied with the QCD sum 
rules [35]- 

In this article, we extend our previous works to study the mass spectrum of the axial- 
vector hidden charmed and hidden bottom tetraquark states in a systematic way with 
the QCD sum rules, and make possible explanation for the nature of the Z + (4430). The 
mass is a fundamental parameter in describing a hadron, whether or not there exist those 
hidden charmed or hidden bottom tetraquark configurations is of great importance itself, 
because it provides a new opportunity for a deeper understanding of the low energy QCD. 
The axial-vector hidden charmed (cc) and hidden bottom (bb) tetraquark states may be 
observed at the LHCb, where the bb pairs will be copiously produced with the cross section 
about 500 fib [36J. 

The hidden charmed and hidden bottom tetraquark states (denoted as Z) have the 
symbolic quark structures: 

Z+ = QQud; Z° = -^QQ(uu - dd); Z~ = QQdu; 
v2 

Zf = QQus; Z~ = QQsu; Z® = QQds; Z s = QQsd; 

Z v = -^=QQ(uu + dd); Z^ = QQss , (1) 
V2 

where the Q denote the heavy quarks c and b. 

We take the diquarks as the basic constituents following Jaffe and Wilczek {371 [35] , and 
construct the axial-vector tetraquark states with the diquark and antidiquark pairs. The 
diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vector 07^75, axial- 
vector C'Jn and tensor Ca^, where C is the charge conjunction matrix. The structures 
C^/fj, and Co^y are symmetric, the structures C75, C and C^^^ are antisymmetric. The 
attractive interactions of one-gluon exchange favor formation of the diquarks in color 
antitriplet 3 C , flavor antitriplet 3f and spin singlet l s [391140], In this article, we assume the 
axial-vector hidden charmed and hidden bottom tetraquark states Z consist of the C75 — 
C^ffi type rather than C — (^7^75 type diquark structures, and construct the interpolating 
currents J M (x) and 77^ (x): 

4 + {x) = e^ k € imn U j(x)C 7 Mx)Q m (x)rCdl(x), 

J Zo( X ) = [ u j(. x ) C l5Qk(x)Qm(xh^Cul(x) - (U -)• d)} , 

r zt {x) = e^ k e mn uj(x)C l5 Q k (x)Q m (x)rCs T n (x), 
J^(x) = e^ k e imn dj(x)C l5 Q k (x)Q m (x)^Csl(x), 

ijk imn 

4 v {x) = —^-[ u j( x )C 75 Q k (x)Q m (x)rCul(x) + ( U ^d)}, 

r z (x) = e^ k e imn sj(x)C 75 Q k (x)Q m (x)^Csl(x), (2) 
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rf z+ (x) = e^e imn U j(x)C^Q k (x)Q m (xh 5 Cd^(x), 

Vzo{x) = -j= [uj (x)C^Q k (x)Q m (x)~f 5 Cul{x) - (u -> d)] , 

rf zt {x) = e^ k e imn uj(x)CYQ k (x)Q m (x) l5 Csl(x), 
Vz o(x) = e^ k e imn dj(x)C^Q k (x)Q m (x) 75 Csl(x), 

AjkAmn 

rfe v (x) = y=- [uj(x)CrQk(x)Q m (x) l5 Cul(x) + (u d)] , 



rt 

J 3 



rfzix) = e^ k e imn sj(x)C^Q k (x)Q m (xh 5 Csl(x), (3) 



where the i, j, k, ■ ■ ■ are color indexes. In the isospin limit, the interpolating currents 
result in three distinct expressions for the spectral densities, which are characterized by 
the number of the s quark they contain. The interpolating currents J^(x) and rj^(x) lead 
to the same expression for the correlation functions IL^ip), for example, 

J z+ ~ ; J z- ~ <- ' J z? ~ <° ; 4+ ~ ^+ ; 

^~^> (4) 

where we use ~ to denote the two interpolating currents lead to the same expression. 
The special superpositions t,J^(x) + (1 — t)r^{x) can't improve the predictions remarkably, 
where t = — 1. In this article, we take the interpolating currents J M (x) for simplicity, 
i.e. t = 1. 

In fact, we can take the colored diquarks as point particles and describe them with the 
scalar S a , pseudoscalar P a , vector V£, axial-vector A® and tensor T® v fields, respectively, 
where the a is the color index, then introduce the SU(3) color interaction. We construct 
the color singlet tetraquark currents with the diquark fields 5°, P a , V£, and T^ u , 
parameterize the nonpertubative effects with the new vacuum condensates (SS), (PP), 
(VV), (AA) and (TT) besides the gluon condensate, and perform the standard procedure 
of the QCD sum rules to study the tetraquark states. The basic parameters such as the 
diquark masses and the new vacuum condensates can be fitted phenomenally. The nonet 
scalar mesons below 1 GeV (the /o(980) and ao(980) especially) are good candidates for the 
tetraquark states, from those tetraquark candidates, we can obtain the basic parameters 
and extend the new sum rules to other tetraquark states. As there are many works to do, 
we prefer another article. 

The article is arranged as follows: we derive the QCD sum rules for the axial-vector 
hidden charmed and hidden bottom tetraquark states Z in Sect. 2; in Sect. 3, we present 
the numerical results and discussions; and Sect. 4 is reserved for our conclusions. 
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2 QCD sum rules for the axial-vector tetraquark states Z 



In the following, we write down the two-point correlation functions H fiu (p) in the QCD 
sum rules, 



IV (p) 



J d 4 xe^(0|T [^(x)jt(O)] |0>, 



(5) 



where the J M (x) denotes the interpolating currents J^ + (x), J^ (x), J^ + (x), etc. 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J^(x) into the correlation functions U^ u (p) to obtain 
the hadronic representation |41|, 142] . After isolating the ground state contribution from 
the pole term of the Z, we get the following result, 



IV (p) 



\2 



p z 



-9/iu + 



+ 



where the pole residue (or coupling) Xz is defined by 



(0\J,(0)\Z(p)) , 



(6) 



(7) 



the denotes the polarization vector. 

After performing the standard procedure of the QCD sum rules, we obtain the following 
six sum rules: 



dspz(s)e 



(8) 



where the Z denote the channels ccqq, ccqs, cess, bbqq, bbqs and bbss respectively; the s° z 
are the corresponding continuum threshold parameters, and the M 2 is the Borel parameter. 
The thresholds Az can be sorted into three sets, we introduce the qq, qs and ss to 
denote the light quark constituents in the axial-vector tetraquark states to simplify the 
notation, A qq - = ^nig, A gS = (2mg + m s ) 2 , A S g = 4(mg + m s ) 2 . The explicit expressions 
of the spectral densities p gq -(s), p q s(s) and p s s(s) are presented in the appendix, where 

l+^l-4mys _ l-^l-4m 2 Q /s _ ^ ( a+/3 ) m 2 

a f ~ 2 ' tti ~ 2 ' ^ ~ Tus^ 



~2 
m Q 



2 > - 2 ) M - as-m^ > ""Q ~ a/3 ' _ a(l-a) ' 

We carry out the operator product expansion to the vacuum condensates adding up to 
dimension-10. In calculation, we take vacuum saturation for the high dimension vacuum 
condensates, they are always factorized to lower condensates with vacuum saturation in 
the QCD sum rules, factorization works well in large A^ c limit. In reality, N c = 3, some 
ambiguities may come from the vacuum saturation assumption. 

We take into account the contributions from the quark condensates, mixed conden- 
sates, and neglect the contributions from the gluon condensate. The gluon condensate 
( a " ) is of higher order in a s , and its contributions are suppressed by very large denom- 
inators comparing with the four quark condensate (qq) 2 (or (ss} 2 ) and would not play any 
significant role, although the gluon condensate ( a " ) has smaller dimension of mass than 
the four quark condensate (qq) 2 (or (ss) 2 ). One can consult the sum rules for the light 
tetraquark states [33j 03], the heavy tetraquark state [31] and the heavy molecular states 
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[45, 46 for example. Furthermore, there are many terms involving the gluon conden- 
sate for the heavy tetraquark states and heavy molecular states in the operator product 
expansion (one can consult Refs.[3lJ US]), we neglect the gluon condensate for simplicity. 

In the special case of the 1"(4660) (as a t///o(980) bound state) and its pseudoscalar 
partner rj' c f o(980), the contributions from the gluon condensate ( a " ) are rather large 
[471 [48] . If we take a simple replacement s(x)s(x) — > (ss) and \u{x)u(x) + d(x)d(x)] —> 
2(qq) in the interpolating currents, the standard heavy quark currents Q(x) / y fM Q(x) and 
Q(x)i^Q(x) are obtained, where the gluon condensate ( a " ) plays an important rule 
in the QCD sum rules The interpolating currents constructed from the diquark- 

antidiquark pairs do not have such feature. There are other interpretations for the y(4660), 
for example, the diquark-antidiquark type charmed baryonium |49| . 

We also neglect the terms proportional to the m u and m^, their contributions are of 
minor importance due to the small values of the u and d quark masses. 

Differentiating the Eq.(8) with respect to ™-, then eliminate the pole residues Xz, we 
can obtain the sum rules for the masses of the Z, 



$l z z ds Pz {s)e F 

3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2)(qq), (qg s aGq) = m 2 (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 0.2) GeV 2 , 
m s = (0.14 ± 0.01) GeV, m c = (1.35 ± 0.10) GeV and m b = (4.8 ± 0.1) GeV at the energy 
scale (i = 1 GeV gH 021 HQ] . 

The Q-quark masses appearing in the perturbative terms are usually taken to be 
the pole masses in the QCD sum rules, while the choice of the tbq in the leading-order 
coefficients of the higher-dimensional terms is arbitrary |5H [52] . The MS mass m c (m 2 ) 

r ( 1\ ~| — 1 

relates with the pole mass rh c through the relation m c (m 2 ) = m c 1 H — + • • ■ 

In this article, we take the approximation m c (m 2 ) ~ rh c without the a s corrections for 
consistency. The value listed in the Particle Data Group is m c (m 2 ) = 1.27^q'ii GeV [53], 
it is reasonable to take rh c = m c (1 GeV 2 ) = (1.35 ± 0.10) GeV. For the b quark, the MS 
mass mb(m 2 ) = 4.201q'q7 GeV ^ e S a P between the energy scale \i = 4.2 GeV and 
1 GeV is rather large, the approximation m& ~ m^m 2 ) ~ m^l GeV 2 ) seems rather crude. 
It would be better to understand the quark masses m c and nib we take at the energy scale 
// 2 = 1 GeV 2 as the effective quark masses (or just the mass parameters). 

In calculation, we also neglect the contributions from the perturbative corrections. 
Those perturbative corrections can be taken into account in the leading logarithmic ap- 
proximations through anomalous dimension factors. After the Borel transform, the effects 
of those corrections are to multiply each term on the operator product expansion side by 

MM2) l2r -< 



the factor, ^Qj! 2 ) " ' w h ere the r j is the anomalous dimension of the interpolat- 

ing current J(x) and the To n is the anomalous dimension of the local operator O n (0). 
We carry out the operator product expansion at a special energy scale /U 2 = 1 GeV 2 , and 

~ 1, such an approximation maybe result in some scale 



set the factor 



a s {M 2 ) 1 2T > 
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dependence and weaken the prediction ability. In this article, we study the axial-vector 
hidden charmed and hidden bottom tetraquark states systemically, the predictions are 
still robust as we take the analogous criteria in those sum rules. 

In the conventional QCD sum rules [41^142] . there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter M 2 and 
threshold parameter sq. We impose the two criteria on the axial- vector heavy tetraquark 
states to choose the Borel parameter M 2 and threshold parameter sq. 

The contributions from the high dimension vacuum condensates in the operator prod- 
uct expansion are shown in Figs. 1-2, where (and thereafter) we use the {qq) to denote the 
quark condensates {qq), {ss) and the {qg s aGq) to denote the mixed condensates {qg s aGq), 
{sg s aGs). From the figures, we can see that the contributions from the high dimension 
condensates are very large and change quickly with variation of the Borel parameter at 
the values M 2 < 2.6 GeV 2 and M 2 < 7.2 GeV 2 in the hidden charmed and hidden bot- 
tom channels respectively, such an unstable behavior cannot lead to stable sum rules, our 
numerical results confirm this conjecture, see Fig. 4. 

At the values M 2 > 2.6 GeV 2 and s > 22 GeV 2 , 23 GeV 2 , 23 GeV 2 , the contributions 
from the {qq) 2 + {qq){qg s &Gq) term are less than 12%, 5%, 2.5% in the channels ccqq, ccqs, 
cess respectively; the contributions from the vacuum condensate of the highest dimension 
{qg s aGq) 2 are less than 2.5%, 2%, 1.5% in the channels ccqq, ccqs, cess respectively; we 
expect the operator product expansion is convergent in the hidden charmed channels. 

At the values M 2 > 7.2 GeV 2 and s > 136 GeV 2 , 138 GeV 2 , 138 GeV 2 , the contri- 
butions from the {qq) 2 + {qq){qg s &Gq) term are less than 10%, 4.5%, 7% in the channels 
bbqq, bbqs, bbss respectively; the contributions from the vacuum condensate of the highest 
dimension {qg s aGq) 2 are less than 5.5%, 4%, 3% in the channels bbqq, bbqs, bbss respec- 
tively; we expect the operator product expansion is convergent in the hidden bottom 
channels. 

In this article, we take the uniform Borel parameter M^ in , i.e. M^ in > 2.6 GeV 2 and 
^min — 7-2 GeV 2 in the hidden charmed and hidden bottom channels respectively. 

In Fig. 3, we show the contributions from the pole terms with variation of the Borel 
parameters and the threshold parameters. The pole contributions are larger than (or 
equal) 46%, 50%, 50% at the value M 2 < 3.2 GeV 2 and s > 22 GeV 2 , 23 GeV 2 , 23 GeV 2 
in the channels ccqq, ccqs, cess respectively, and larger than (or equal) 48%, 50%, 50% 
at the value M 2 < 8.2 GeV 2 and s > 136 GeV 2 , 138 GeV 2 , 138 GeV 2 in the channels 
bbqq, bbqs, bbss respectively. Again we take the uniform Borel parameter M^ aa ,, i.e. 
M^ ai < 3.2 GeV 2 and M^ aa . < 8.2 GeV 2 in the hidden charmed and hidden bottom 
channels respectively. 

In this article, the threshold parameters are taken as sq = (23±l)GeV 2 , (24±l)GeV 2 , 
(24 ± 1) GeV 2 , (138 ± 2) GeV 2 , (140 ± 2) GeV 2 , (140 ± 2) GeV 2 in the channels ccqq, ccqs, 
cess, bbqq, bbqs, bbss respectively; the Borel parameters are taken as M 2 = (2.6-3.2) GeV 2 
and (7.2 — 8.2) GeV 2 in the hidden charmed and hidden bottom channels respectively. In 
those regions, the pole contributions are about (46 — 74)%, (50 — 77)%, (50 — 77)%, 
(48 — 67)%, (50 — 69)%, (50 — 69)% in the channels ccqq, ccqs, cess, bbqq, bbqs, bbss 
respectively; the two criteria of the QCD sum rules are fully satisfied |41l 142] . 

From Fig.3, we can see that the Borel windows M^ ax — M^ in change with variations of 
the threshold parameters sq. In this article, the Borel windows are taken as 0.6 GeV 2 and 
1.0 GeV 2 in the hidden charmed and hidden bottom channels respectively; they are small 
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enough. If we take larger threshold parameters, the Borel windows are larger and the 
resulting masses are larger, see Fig. 4. In this article, we intend to calculate the possibly 
lowest masses which are supposed to be the ground state masses by imposing the two 
criteria of the QCD sum rules. 

If we take analogous pole contributions, the interpolating current with more s quarks 
requires slightly larger threshold parameter due to the SU(3) breaking effects, see Fig. 3. 
In the channels QQqs and QQss, the SU(3) breaking effects on the threshold parameters 
are tiny, we take uniform threshold parameters in those channels. In Fig. 4, we plot the 
axial- vector tetraquark state masses Mz with variation of the Borel parameters and the 
threshold parameters. Naively, we expect the tetraquark state with more s quarks will have 
larger mass. In calculations, we observe that the possibly lowest masses of the axial- vector 
heavy tetraquark states QQqs and QQss are almost the same. 

Taking into account all uncertainties of the relevant parameters, finally we obtain the 
values of the masses and pole resides of the axial-vector tetraquark states Z, which are 
shown in Figs. 5-6 and Table 1. In Table 1, we also present the masses of the scalar hidden 
charmed and hidden bottom tetraquark states obtained in our previous works |32|, 133] . 

In this article, we calculate the uncertainties 5 with the formula 



i ^ ' 

where the / denote the hadron mass Mz and the pole residue \z , the x% denote the relevant 
parameters m c , m^, {qq), {ss), ■ ■ ■ . As the partial derivatives Jj- are difficult to carry out 

analytically, we take the approximation (xi — Xi) 2 ~ [f(xi ± Axi) — f(x~i)] 2 in the 

numerical calculations. 

From Table 1, we can see that the uncertainties of the masses Mz are rather small 
(about 4% in the hidden charmed channels and 2% in the hidden bottom channels), while 
the uncertainties of the pole residues Xz are rather large (about 20%). The uncertainties 
of the input parameters ({qq), {ss), {sg s aGs), {qg s aGq), m s , m c and nib) vary in the range 
(2 — 25)%, the uncertainties of the pole residues Xz are reasonable. We obtain the squared 
masses M| through a fraction, the uncertainties in the numerator and denominator which 
originate from a given input parameter (for example, {ss), {sg s aGs)) cancel out with each 
other, and result in small net uncertainty. 

The SU (3) breaking effects for the masses of the axial-vector hidden charmed and hid- 
den bottom tetraquark states are buried in the uncertainties. Naively, we expect the axial- 
vector and vector diquarks have larger masses than the corresponding scalar diquarks, 
and the masses of the tetraquark states have the hierarchy Mq^^-c^ > ^C^-C^ > 
Mc 75 _c* 75 , because the attractive interactions of one-gluon exchange favor formation of 
the diquarks in color antitriplet 3 C , flavor antitriplet 3/ and spin singlet l s [39} 140] . From 
Table 1, we can see that it is not the case. 

In the conventional QCD sum rules, we usually consult the experimental data in choos- 
ing the Borel parameter M 2 and the threshold parameter sq. If the mass spectrum of the 
axial-vector tetraquark states are well known, we can denote the ground state, the first 
excited state, the second excited state, the third excited state, . . . , as Z, Z', Z", Z'" , 
. . . . The critical thresholds for emergence of those excited tetraquark states are Tz>, 
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Tz" , Tz"i , . . . , respectively. The threshold parameter so should take values in the region 
(Mz + Tz) 2 < so < T Z '. However, the present experimental knowledge about the phe- 
nomenological hadronic spectral densities of the multiquark states is rather vague, even 
the existence of the multiquark states is not confirmed with confidence, and no knowledge 
about either there are high resonances or not. 

Taking into account the two criteria (pole dominance and convergence of the operator 
product expansion) of the QCD sum rules, we can obtain the possibly lowest threshold 
parameter so, which is denoted as s^ in . I* 1 this article, we take the value s^ n and make 
crude estimations for the ground state masses. 

The values of the s^ in in different channels maybe smaller (or larger) than T Z ', or 
even smaller than (Mz + Tz) 2 , the two criteria of the QCD sum rules alone cannot always 
warrant satisfactory threshold parameters and Borel windows. For example, the nonet 
scalar mesons below 1 GeV (the /o(980) and ao(980) especially) are good candidates for 
the tetraquark states |38|, [541 [55] . The two criteria of the QCD sum rules result in the 
threshold parameters so 3> (Mj 0/ / ao + Tf / ao ) 2 , the contributions of the excited states are 
already included in if there are any, and we have to resort to "multi-pole + continuum 
states" to approximate the phenomenological spectral densities. If we insist on the "one- 
pole + continuum states" ansatz, no reasonable Borel window can be obtained, although it 
is not an indication non-existence of the light tetraquark states (For detailed discussions 
about this subject, one can consult Refs. [3Tl I56j). The QCD sum rules is just a QCD 
model. 

In the channel ccqq, the threshold parameter s^ in leads to the mass M C c q q = (4.32 ± 

0. 18) GeV, which is consistent with the experimental data Mz = (4433 ± 4 ± 2) MeV or 
4443^12^13 MeV from the Belle collaboration within uncertainty [TUB]- The experimental 
value is (M z + T z ) 2 < 22.5 GeV 2 , the lower bound of the s° min = (22 - 24) GeV 2 is 
smaller than (Mz + T^) 2 , we have to postpone the s^ in to larger values. If we take 
so = (26 ± 1) GeV 2 , the prediction Mz = (4.44 ± 0.19) GeV is in excellent agreement with 
experimental data, see Fig. 7. In Fig. 8, we present the corresponding pole residue, from 
the figure, we can see that larger threshold parameter result in larger pole residue. 

The predictions of the QCD sum rules favor the scenario of the Z + (4430) as an axial- 
vector tetraquark state, the Z + (4430) can be tentatively identified as an axial- vector 
tetraquark state. In other channels, the heavy axial-vector tetraquark states exist in 
nature maybe have larger masses than those theoretical predictions presented in Table 

1. On the other hand, the upper bound of the threshold parameter so = (25 — 27) GeV 2 
maybe larger than the critical threshold Tz>, so we identify the Z + (4430) as an axial- vector 
tetraquark state tentatively, not confidently. 

In Refs.[3T| l32"l [33], we observe that the meson Z(4250) may be a scalar tetraquark 
state (ccud), irrespective of the C^^ — C^ type and the C75 — C75 type diquark structures, 
the decay Z(4250) — > ir + Xd can take place with the Okubo-Zweig-Iizuka super-allowed 
"fall-apart" mechanism, which can take into account the large total width naturally. In 
the present case, the decay Z(4430) — > ip'tr can also take place with the Okubo-Zweig- 
Iizuka super-allowed "fall-apart" mechanism, which can take into account the large total 
width (T z = 45±^t^MeV or 107±||±^| MeV) naturally [3[8]. 

In this article, we take the simple pole + continuum approximation for the phenomeno- 
logical spectral densities. In fact, such a simple approximation has shortcomings. In the 
case of the non-relativistic harmonic-oscillator potential model, the spectrum of the bound 
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states (the masses E n and the wave functions ^> n (x)) and the exact correlation functions 
(and hence its operator product expansion to any order) are known precisely. The non- 
relativistic harmonic-oscillator potential ^mw¥ is highly non-perturbative, one suppose 
the full Green function satisfies the Lippmann-Schwinger operator equation and may be 
solved perturbatively. We can introduce the Borel parameter dependent effective threshold 

parameter z e f f (M) = ui [zq + z\ + ^2 H ] and fit the coefficients to reproduce 

both the ground energy Eq and the pole residue Rq = Vl/p (0)^0(0), the phenomenological 
spectrum density can be described by the perturbative contributions well above the effec- 
tive continuum threshold z e ff(M), or reproduce the ground energy Eq only and take the 
pole residue R as a calculated parameter, there exists a solution for the effective contin- 
uum threshold z e fj(M) which precisely reproduces the exact ground energy Eq for any 
value of the pole residue R within the range R = (0.7 — 1.15)12q in the limited fiducial 
Borel window, the value of the pole residue R extracted from the sum rule is determined 
to a great extent by the contribution of the hadron continuum [57j. There maybe systemic 
uncertainties out of control. 

In the real QCD world, the hadronic spectral densities are not known exactly. In the 
present case, the ground states have not been observed yet, except for the possible axial- 
vector tetraquark state candidate Z + (4430). So we have no confidence to introduce the 
Borel parameter dependent effective threshold parameter s e ff(M) = sq+sijj?+S2 • • 
and approximate the phenomenological spectral densities with the perturbative contribu- 
tions above the effective continuum threshold s e ff(M) accurately. Furthermore, the pole 
residues (or the couplings of the interpolating currents to the ground state tetraquark) Xz 
are not experimentally measurable quantities, and should be calculated by some theoret- 
ical approaches, the true values are difficult to obtain, which are distinguished from the 
decay constants of the pseudoscalar mesons and the vector mesons, the decay constants 
can be measured with great precision in the leptonic decays (in some channels). 

The spectrum of the bound states in the non-relativistic harmonic-oscillator potential 
model are of the Dirac 5 function type, we can choose z e ff < E\, while in the case of 
the QCD, the situation is rather complex, the effective continuum thresholds s e tf(M) 
maybe overlap with the first radial excited states, which are usually broad. For example, 
in the pseudoscalar channels, the widths of the tt, 7r(1300), 7r(1800), ••• are ~ OGeV, 
(0.2-0.6) GeV, (0.208±0.012) GeV, • • • respectively, while the widths of the K, #(1460), 
#(1830), ••• are ~ OGeV, ~ (0.25 - 0.26) GeV, ~ 0.25GeV, ••• respectively pj. In 
this article, we prefer (or have to choose) the simple pole + continuum approximation, 
and cannot estimate the unknown systemic uncertainties of the QCD sum rules before 
the spectral densities in both the QCD and phenomenological sides are known with great 
accuracy. 

In Ref. |58j . Lucha, Melikhov and Simula use the correlation function of the pseu- 
doscalar current J§(x) = (mb + m u )q(x)i^b(x) to illustrate a Borel-parameter-dependent 
effective continuum threshold can reduce considerably the (unphysical) dependence of the 
extracted bound-state mass and the decay constant on the Borel parameter. In the present 
case, we have no experimental data for the masses and pole residues of the tetraquark 
states to take as a guide and apply the x 2 minimization by adjusting the effective thresh- 
old parameters. On the other hand, the Borel-parameter-dependent effective continuum 
thresholds maybe overlap with the (Mz + ^z) 2 , Tz>, Tz», we prefer the Borel- 
parameter-independent threshold parameters, although the Borel-parameter-dependent 



9 



effective threshold parameters maybe smear some dependence on the Borel parameter. 
From Figs. 5-8, we can see that the dependence of the masses and pole residues on the 
Borel parameters in the Borel windows are rather mild. 

The central values of our predictions are much larger than the corresponding ones 
from a relativistic quark model based on a quasipotential approach in QCD [59j[60]. In 
Refs.[59, 60], Ebert et al take the diquarks as bound states of the light and heavy quarks in 
the color antitriplet channel, and calculate their mass spectrum using a Schrodinger type 
equation, then take the masses of the diquarks as the basic input parameters, and study the 
mass spectrum of the heavy tetraquark states as bound states of the diquark-antidiquark 
system. In Refs.[6l] [62l [63], Maiani et al take the diquarks as the basic constituents, 
examine the rich spectrum of the diquark-antidiquark states from the constituent diquark 
masses and the spin-spin interactions, and try to accommodate some of the newly ob- 
served charmonium-like resonances not fitting a pure cc assignment; furthermore, the 
corresponding bottom tetraquark states are also studied with the same method [53]. The 
predictions depend heavily on the assumption that the light scalar mesons ao(980) and 
/o(980) are tetraquark states, the basic parameters (constituent diquark masses) are es- 
timated thereafter. In the conventional quark models, the constituent quark masses are 
taken as the basic input parameters, and fitted to reproduce the mass spectra of the 
conventional mesons and baryons. However, the present experimental knowledge about 
the phenomenological hadronic spectral densities of the tetraquark states is rather vague, 
whether or not there exist tetraquark states is not confirmed with confidence, and no 
knowledge about the high resonances. The predicted constituent diquark masses can not 
be confronted with the experimental data. 

In Refs.[2T | l6"0" l l62]. the X(3872) and Z(4430) are taken as the ground state axial-vector 
and first radially excited axial- vector tetraquark states respectively. In Ref. [65j . Matheus 
et al study the X(3872) as a tetraquark state with J PC = 1 ++ using QCD spectral sum 
rules, the prediction is consistent with the experimental data within uncertainty. The dis- 
crepancy between the predictions of Ref. [65] and the present work (analogous interpolating 
currents are chosen in those works) mainly originates from the high dimensional vacuum 
condensates (qg s aGq) 2 which are neglected in Ref. [65]. The condensates {qg s aGq) 2 are 
counted as O(jfe), and the corresponding contributions are greatly enhanced at small M 2 , 
and result in rather bad convergent behavior in the operator product expansion, we have to 
choose larger Borel parameter M 2 . We insist on taking into account the high dimensional 
vacuum condensates, as the interpolating current consists of a light quark- antiquark pair 
and a heavy quark-antiquark pair, one of the highest dimensional vacuum condensates is 
(qq) 2 x ), we have to take into account the condensates (qg s aGq) 2 for consistence. 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. In 
proton-proton collisions at yfs = 14TeV, the bb cross section is expected to be ~ 500/x6 
producing 10 12 bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 x 10 32 cm _2 sec _1 [36]. The axial- vector tetraquark states predicted in the present 
work may be observed at the LHCb, if they exist indeed. We can search for the axial- 
vector hidden charm tetraquark states in the DD*, DD*, D S D*, D S D*, J/ipir, J/rpK, 
J/iprj, V 7 ' 71 "! i^K-i ' ' ' invariant mass distributions and search for the axial-vector hidden 
bottom tetraquark states in the BB*, BB*, B S B*, B S B*, Tvr, TK, Trj, Vir, T'K, T% 
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Figure 1: The contributions from different terms with variation of the Borel parameter 
M 2 in the operator product expansion. The (I) and (II) denote the contributions from 
the (qq) 2 + {qq}(qg s ^Gq} term and the (qg s aGq) 2 term respectively. The A, B and C 
denote the channels ccqq, ccqs and cess respectively. The notations a, [3, 7, A, r and 
p correspond to the threshold parameters s = 21 GeV 2 , 22 GeV 2 , 23 GeV 2 , 24 GeV 2 , 
25 GeV 2 and 26 GeV 2 respectively. 
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Figure 2: The contributions from different terms with variation of the Borel parameter 
M 2 in the operator product expansion. The (I) and (II) denote the contributions from 
the (qq) 2 + {qq){qg s &Gq) term and the (qg s aGq) 2 term respectively. The A, B and C 
denote the channels bbqq, bbqs and bbss respectively. The notations a, f3, 7, A, r and p 
correspond to the threshold parameters so = 132 GeV 2 , 134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 
140 GeV 2 and 142 GeV 2 respectively. 
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Figure 3: The contributions of the pole terms with variation of the Borel parameter 
M 2 . The A, B, C, D, E and F denote the channels ccqq, ccqs, cess, bbqq, bbqs and bbss 
respectively. The notations a, [3, 7, A, r and p correspond to the threshold parameters so = 
21 GeV 2 , 22 GeV 2 , 23 GeV 2 , 24 GeV 2 , 25 GeV 2 and 26 GeV 2 respectively in the hidden 
charmed channels; while they correspond to the threshold parameters sq = 132 GeV 2 , 
134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 140 GeV 2 and 142 GeV 2 respectively in the hidden bottom 
channels. 
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Figure 4: The masses of the axial-vector tetraquark states with variation of the Borel 
parameter M 2 . The A, B, C, D, E and F denote the channels ccqq, ccqs, cess, bbqq, 
bbqs and bbss respectively. The notations a, (3, 7, A, r and p correspond to the threshold 
parameters s = 21GeV 2 , 22GeV 2 , 23GeV 2 , 24GeV 2 , 25GeV 2 and 26GeV 2 respectively 
in the hidden charmed channels; while they correspond to the threshold parameters sq = 
132 GeV 2 , 134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 140 GeV 2 and 142 GeV 2 respectively in the 
hidden bottom channels. 
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Figure 5: The masses of the axial- vector tetraquark states with variation of the Borel 
parameter M 2 . The A, B, C, D, E and F denote the channels ccqq, ccqs, cess, bbqq, bbqs 
and bbss respectively. 
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Figure 6: The pole residues of the axial-vector tetraquark states with variation of the 
Borel parameter M 2 . The A, B, C, D, E and F denote the channels ccqq, ccqs, cess, 
bbqq, bbqs and bbss respectively. 
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Figure 8: The pole residue of the Z(4430) with variation of the Borel parameter M 2 . 
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tetraquark states 


c 75 - c ltx 




C% - 


C75 - c 75 


ccqq 


4.32 ±0.18 


2.92 ± 0.60 


4.36 ±0.18 


4.37 ±0.18 


ccqs 


4.41 ±0.16 


3.51 ±0.70 




4.39 ±0.16 


cess 


4.40 ±0.16 


3.51 ±0.70 


4.45 ±0.16 


4.44 ±0.16 


bbqq 


11.27 ±0.20 


1.48 ± 0.34 


11.14 ±0.19 


11.27 ±0.20 


bbqs 


11.38 ±0.18 


1.80 ±0.43 




11.33 ±0.16 


bbss 


11.34 ±0.16 


1.78 ±0.41 


11.23 ±0.16 


11.31 ±0.16 


Z(4430) 


4.44 ±0.19 


3.94 ± 0.71 







Table 1: The masses and the pole residues of the axial- vector tetraquark states. The 
masses are in unit of GeV and the pole residues are in unit of 10 -2 GeV 5 and 10 _1 GeV 5 
in the channels cc and bb respectively. 

• • • invariant mass distributions. 

4 Conclusion 

In this article, we construct the scalar-diquark-axial-vector-antidiquark type currents to 
interpolate the axial-vector tetraquark states, and study the mass spectrum of the axial- 
vector hidden charmed and hidden bottom tetraquark states with the QCD sum rules in a 
systematic way. In calculations, we take into account the contributions from the vacuum 
condensates adding up to dimension 10 in the operator product expansion, and neglect 
the gluon condensates as their contributions are supposed to be very small. The mass 
spectrum are calculated by imposing the two criteria (pole dominance and convergence of 
the operator product expansion) of the QCD sum rules. Our numerical result M z ^ 430 ^ = 
(4.44 ± 0.19) GeV is in excellent agreement with the experimental data Mz = (4433 ± 
4 ± 2)MeV or 4443^^^13 Me ^ ^ rom tne Bene collaboration, which indicates that the 
Z + (4430) can be identified as the axial- vector tetraquark state tentatively. Considering the 
light-flavor SU (3) symmetry and the heavy quark symmetry, we make predictions for other 
axial-vector tetraquark states, the predictions can be confronted with the experimental 
data in the future at the LHCb or the Fermi-lab Tevatron. 

Appendix 

The spectral densities p qq -(s), p q s(s) and p S s(s) at the level of the quark-gluon degrees of 
freedom: 
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